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Objectives 


Concepts  of  Space-Time  Coding  (STC)  for  communication  signals  have  been  well  researched  and  implemented  in 
the  domain  of  MIMO  telecommunication  systems.  The  aim  of  this  project  has  been  to  explore  how  far  such  signal 
design  aspects  could  be  carried  over  to  MIMO  radar  systems  for  exploiting  the  additional  degrees  of  freedom 
introduced  by  STC.  Examples  of  MIMO  radars  are  phased  arrays  and  netted  radars  in  the  context  of  network  centric 
warfare.  The  fundamental  problem  here  is  to  answer  what  is  an  optimal  radar  receiver  in  the  context  of  signal 
processing  and  how  to  solve  the  mathematical  intricacies  of  coding.  Once  this  issue  is  addressed,  one  would  like  to 
examine  the  performance  of  such  a  receiver  for  various  radar  related  applications,  for  example  for  target  tracking. 

Status  of  effort 

We  have  presented  in  the  Conferences  all  our  research  results.  We  will  be  extending  this  work  in  further  depth  and 
have  added  a  section  in  this  regard  on  future  direction  at  the  end  of  this  report  for  perusal.  It  is  our  objective  to 
consolidate  all  existing  work  in  an  archival  journal  paper  to  be  submitted  in  2006.  There  remain  several  remaining 
unsolved  problems,  mainly  in  respect  to  the  solution  to  the  fundamental  optimisation  problems  arising  from  the 
work. 

Abstract 

The  conceptual  framework  for  designing  transmit  signals  for  MIMO  radar  systems  in  terms  of  the  techniques  of  STC 
as  developed  for  communication  signals  has  been  formulated.  This  is  achieved  by  constructing  a  generalized 
matched  filter  in  analogy  of  the  usual  case.  We  have  studied  in  detail  the  corresponding  maximum  likelihood 
receiver  and  shown  this  to  be  optimal  in  accordance  with  the  Cramer-Rao  lower  bound  criterion.  We  have  developed 
the  mathematics  for  defining  and  generating  the  coding  coefficients.  Tracking  of  a  single  target  has  been  treated 
within  this  framework  in  terms  of  state  equations.  A  tracker  defined  by  Kalman  predictor  has  been  introduced.  A 
channel/target  model  involving  both  linear  and  nonlinear  dynamics  has  been  studied.  We  have  briefly  touched  on  the 
implication  of  STC  for  transmit  beam  forming. 
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Interactions 


As  noted  above  the  results  of  the  research  have  been  presented  in  the  First  and  Second  International  Conference  on 
Waveform  Diversity  and  Design  (Dr.  P.S.Ray)  and  in  the  Asilomar  Conference  (Prof.  L.B. White)  which  was  an 
invited  session.  There  has  been  considerable  scientific  interaction  with  international  experts  at  the  Conferences  as 
described  in  the  following. 

First,  Dr.  Ray  had  a  brief  discussion  with  Dr.  Mike  Wicks  ( AFRL)  who  advised  to  keep  him  up-to-date  with  our 
further  progress  in  the  research.  Fie  has  assured  his  support  for  any  purpose  when  required.  He  endorsed  a  “short” 
visit  by  Dr.  Ray  to  Hanscom  (AFRL)  in  order  to  create  an  opportunity  for  his  participation  at  the  upcoming  SAM 
2006  Workshop  (Sensor  Array  and  Multi-Channel  Signal  Processing).  This  is  subject  to  the  immediate  availability 
of  funding  from  AFRL,  which  is  being  explored  by  Dr.  Rangaswamy  (AFRL). 

Dr.  Marshall  Greenspan,  Technical  Director  of  Northrop  Grumman  Corporation,  indicated  to  Dr.  Ray  that  he 
intended  to  draw  the  attention  of  some  group  in  the  Ohio  State  University  to  the  expertise  in  Adelaide  in  MIMO 
coding  techniques  for  possible  future  collaboration.  They  are  planning  some  activity  on  MIMO  systems  for  short 
range  applications. 

A  promising  development  has  come  from  Professor  Chris  Baker  in  University  College  London.  He  wanted  us  to 
submit  a  program  of  research  on  coding  of  MIMO  radars  for  collaborative  work.  His  group  will  undertake 
experimental  investigations  in  this  area  as  they  have  hardware  designing  expertise.  We  are  following  up  this  in  due 
course. 

Following  the  Asilomar  conference  in  2004,  Dr  Rick  Blum  from  Lehigh  University,  USA  invited  Prof  White  to  be  an 
associate  investigator  on  a  MURI  application  for  MIMO  radar  system  design.  This  application  was  unsuccessful  but 
did  serve  the  purpose  of  familiarising  the  research  teams  at  each  institution  about  each  others'  work. 

New 

The  work  undertaken  here  falls  under  the  category  of  “enabling  research”  and  as  such  it  is  early  to  assess  at  present 
its  impact  in  the  field.  We  can  only  vouch  for  its  originality  by  our  publications.  We  believe  that  our  work  is  the  first 
to  use  a  model  based  approach  to  extend  earlier  work  by  David  Kershaw  and  Rob  Evans  (IEEE  Trans,,  on 
Information  Theory,  1994)  to  the  MIMO  case.  We  also  believe  that  we  are  the  first  to  link  the  use  of  unitary 
transmitter  codes  to  obtain  a  simple  correlation  based  receiver  which  achieves  close  to  maximum  likelihood 
performance.  This  result  has  the  added  implication  that  obtaining  transmit  directionality  by  use  of  non-unitary  codes 
may  lead  to  either  the  requirement  for  an  excessively  complicated  optimal,  or  perhaps  poorly  performing  sub-optimal 
MIMO  receiver. 

Honors/Awards 

This  research  program  is  at  its  infancy  at  the  present  stage  and  as  such  this  criterion  seems  not  to  be  applicable. 

Archival  Documentation 

The  research  publications  in  refereed  Conferences  are  enclosed. 

Software  and/or  Hardware 

We  have  developed  experimental  software  in  matlab  to  simulate  a  single  target  MIMO  transmitter  and  receiver 
baseband  system. 


Future  Directions 


At  this  stage  I  will  note  the  following  areas  for  immediate  investigation: 

We  need  to  sharpen  our  results  for  the  tracking  application.  A  preliminary  study  was  undertaken  in  this 
direction  in  the  paper  presented  at  the  Edinburgh  Conference  (2004).  In  particular,  there  are  issue  of  the 
effect  of  Doppler  side-lobes  on  the  receiver  parameter  estimation  which  we  need  to  address  properly.  Also, 
it  is  not  clear  whether  the  resulting  code  optimisation  problem  eeven  has  a  feasible  solution,  and  some 
imposed  code  constraints  may  need  to  be  relaxed. 

We  have  just  touched  in  our  present  work  the  implication  of  MIMO  coding  for  the  transmit  beam  forming. 
We  need  to  develop  this  aspect  fully.  The  incorporated  explicit  transmitter  beamforming  into  our  system 
model  has  the  advantage  of  permitting  transmitter  bean  control  whilst  still  retaining  near  ML  receiver 
performance  enabled  by  the  use  of  orthogonal  transmitter  codes. 

The  questions  raised  in  the  conferences  prompt  us  to  take  into  account  clutter  analysis  for  the  receiver 
performance,  a  significant  issue  in  the  radar  context.  We  have  neglected  this  so-far  in  our  work,  as  we  have 
focussed  on  airborne  targets. 

We  have  as  yet  to  deal  with  multiple  target  tracking,  where  we  conjecture  that  the  additional  degrees  of 
freedom  offered  by  MIMO  systems  will  lead  to  considerable  benefits. 


The  fundamentals  of  this  research  fall  in  the  category  of  statistical  signal  processing.  We  have  already  initiated 
some  discussions  with  Professor  Eric  Moulines  (ENST,  Paris)  in  this  context.  We  intend  to  introduce  Monte  Carlo 
techniques  for  multiple  target  simulations  and  the  associated  optimization  technique  in  future.  Profesor  Moulins’s 
group  has  international  expertise  in  this  area.  Prof  White  will  be  hosted  by  Prof  Moulines  during  his  sabbatical  in  the 
first  half  of  2007.  We  will  be  seeking  funding  to  permit  Dr  Ray  to  visit  for  an  extended  period  so  that  work  on 
multiple  target  tracking  can  be  advanced. 
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Abstract —  The  paper  is  the  first  in  a  series  of  three  which 
address  the  problem  of  radar  transmitter  and  receiver  de¬ 
sign  for  MIMO  (multiple  transmit  and  receive  antennae) 
tracking  systems.  This  paper  considers  the  receiver  design 
problem.  The  method  described  consists  of  a  generalised 
matched  filter  receiver  which  produces  approximate  maxi¬ 
mum  likelihood  estimates  of  the  target’s  delay  and  doppler, 
and  the  spatial  wavenumber  (angle  of  arrival)  and  its  time 
derivative.  The  paper  argues  that  this  approach  is  novel 
in  the  sense  that  it  appropriately  combines  the  usual  single 
antenna  matched  filter  receiver,  and  the  beamforming  re¬ 
quired  for  the  multiple  antenna  case.  Simulation  results  are 
used  to  assess  the  performance  of  the  proposed  method. 

I.  Introduction 

There  has  been  considerable  recent  interest  in  the  use  of 
transmit  and  receive  antenna  arrays  to  improve  the  detec¬ 
tion  and  tracking  performance  of  radar  systems,  in  partic¬ 
ular  the  recent  emergence  of  a  number  of  dedicated  work¬ 
shops  addressing  this  problem.  Whilst  the  waveform  de¬ 
sign  problem  has  been  considered  for  single  transmit  an¬ 
tenna  tracking  systems  (see  eg  [2]),  there  has  been  little 
attention  paid  to  the  corresponding  problem  for  MIMO 
(Multiple-Input,  Multiple-Output)  systems.  Recently,  the 
corresponding  problem  for  MIMO  communications  systems 
has  been  considered  [3],  although  there  are  significant  dif¬ 
ferences  between  the  communications  and  radar  problems. 

The  use  of  antenna  arrays  permits  not  only  the  estima¬ 
tion  of  the  targets’  range  and  radial  velocity,  but  also  its 
spatial  wavenumber  which  is  related  to  the  target  position. 
The  time  derivative  of  the  wavenumber  can  also  be  esti¬ 
mated  giving  information  on  the  angular  velocity  of  the 
target  with  respect  to  the  antenna  array.  An  approxima¬ 
tion  to  the  MIMO  antenna  data  likelihood  for  a  single  tar¬ 
get  model,  leads  to  a  generalised  matched  filter  receiver 
architecture  consisting  of  N  matched  filters  for  each  of  the 
N  receiver  antennae.  The  outputs  of  these  N'2  matched 
filters  are  coherently  combined  to  yield  a  statistic  which 
can  be  maximised  to  yield  the  required  estimates.  Due  to 
the  approximation  made  in  deriving  the  statistic,  bias  will 
generally  be  present  in  the  estimates,  especially  for  small 
antenna  arrays.  Local  maxima  can  also  be  a  problem,  es¬ 
pecially  when  estimating  the  velocity  parameters. 

In  this  paper,  we  have  made  no  attempt  to  choose  the 
transmitted  signal  waveforms  in  order  to  optimise  the  re¬ 
ceiver  performance.  In  [6],  we  derive  approximate  Cramer- 
Rao  bounds  for  the  receiver  presented  here,  and  show  how 
the  transmitted  waveforms  can  be  chosen  to  minimise  a 


weighted  mean-square  error  measure  for  the  estimated  pa¬ 
rameters.  In  [7],  the  tracking  problem  is  addressed,  includ¬ 
ing  adaptive  waveform  design  to  optimise  tracking  perfor¬ 
mance. 

II.  Signal  Model 

Consider  an  linear  equispaced  array  of  N  sensors  located 
along  the  y  axis  as  shown  in  figure  1.  The  position  of  the 
n-th  sensor  is  (n  —  1  )dy  where  d  is  the  element  spacing  and 
y  is  the  unit  vector  in  the  y  direction.  We  will  transmit  suc¬ 
cessive  blocks  of  M  pulses,  each  of  duration  T  seconds  from 
each  antenna.  These  signals  have  the  form  (after  upmixing 
on  carrier  of  frequency  c o  (rad/sec)) 

M-l 

Sn,fc(<)  =  elwt  22  Xn,m(k)  gn{t  -  ( kM  +  m)T)  ,  (1) 

m= 0 

for  n  =  0, . . .  ,iV  —  1,  where  gn  is  the  unit  energy  pulse 
shaping  function  for  transmitter  n  supported  on  [0,T]  and 
Xn,m(k)  is  a  complex  code  sequence  transmitted  by  trans¬ 
mitter  n  for  signal  block  k.  The  baseband  quadrature  re¬ 
ceived  signal  at  sensor  n  due  to  a  single  target  reflection  is 
given  by 

N- 1 

rn,k{t)  =  a(k)  (t  -  (■ rn(k )  +  n(k))) 

e=o 

e-Mrn(k)+rdk)) +  Wn^  (2) 

where  u(k)  is  an  unknown  complex  quantity  incorporating 
the  effects  of  medium  attenuation  and  target  reflectivity. 
These  variables  are  modelled  as  zero  mean  complex  Gaus¬ 
sian  process  independent  from  block  to  block.  Here  r„(fc) 
and  vn{k)  are  the  (one-way)  delay  and  doppler  from  sensor 
n  to  the  target  over  the  signal  block  k,  and  wn^{t)  is  com¬ 
plex  Gaussian  zero  mean  white  receiver  noise.  We  assume 
that  the  duration  of  each  signal  block  is  small  compared 
to  the  target  motion  is  sufficiently  slow  with  respect  to  the 
block  length  so  that  the  delay  and  doppler  are  constant 
over  each  signal  block. 

We  assume  that  there  is  a  single  target  located  at  polar 
co-ordinates  (ro(fc),  0o(k),  (j)o{k))  during  signal  block  k  with 
respect  to  the  reference  sensor  located  at  the  origin.  The 
delay  r„(fc)  and  doppler  shift  vn{k)  as  seen  by  sensor  n  are 
related  for  each  n  to  the  reference  delay  and  doppler  seen 
by  the  sensor  at  the  origin  (to  linear  terms,  and  assum¬ 
ing  that  the  target  distance  is  much  larger  than  the  array 
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Fig.  1.  Array  and  Target  Geometry 


We  now  describe  the  form  of  the  generalised  matched 
filter.  The  log  likelihood  function  associated  with  the  mea¬ 
surements  from  the  N  receive  sensors  for  block  k  is 
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.  N-l 

E 

n— 0 


N—l 

rn,k(t )  sjtk  (t 

3=0 


Tn{k) 


Tj(k)) 


e~iu}{Tn(k)+Tt(k))  eit{vn{k)  +  Vj(k)) 


(St) 


where  Tn(k)  and  vn(k)  are  the  (one-way)  delay  and  doppler 
associated  with  sensor  n.  Such  a  function  is  difficult  to 
maximise  over  the  unknown  parameters  due  to  the  cross 
terms  present  in  this  superimposed  signal  model.  Thus  we 
neglect  these  cross  terms  and  use  instead  the  statistic  1 
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N-l 

=  ^  a  Re 

n,j=0 

X  e-Mr„(fc)  +  rf(fc))  e-it(un(k)+uj(k))  _^N/2  a2  Ef_  (6) 


rn,k(t) 


sj,k(t~T^(k)  ~TAk)) 


dimension)  by 

r -„(*)  = 

c 

undp0(k) 

vn(k)  =  u0(k) - , 

c 

where 


To(k) 

V0(k) 


r0(k) 

c 

ur0(k) 

c 


(3) 


(4) 


are  the  one-way  delay  and  doppler  to  the  reference  sen¬ 
sor.  The  wavenumber  variable  for  this  array  geometry  is 
Po(k)  =  sin(6)o(fc))  sin (cj)0(k)).  The  target  state  variables 
are  denoted  9(k)  =  [r0(fc),  vp(k),  po(k),  Po(k)\- 

III.  Receiver  Architecture 
The  overall  system  architecture  is  shown  in  figure  2. 


Waveform  code  for  next  block  of  pulses 


Fig.  2.  Radar  Tracker  System  Architecture 


In  this  paper,  we  only  address  the  design  of  the  gener¬ 
alised  matched  filter,  with  performance  evaluation  (wrt  the 
Cramer- Rao  bound)  and  code  design  discussed  in  [6] .  The 
code  design  problem  for  tracking  systems  is  addressed  in 

[7]- 


where  Ek  is  the  total  transmitted  energy  of  the  signal  dur¬ 
ing  block  k,  and  a  =  a  e'A  Maximising  directly  over  a  and 
4>  yields  the  (approximate)  MLEs 


a 


where 


1 

~NE~k 


£*!>> 


arctan 


_ReE  n,jXnj(0)\ 


(7) 


Xnj(0)  =  J  rn,k(t)  s*>k  (t  -  Tn(k)  -  Tj(k)) 

x  e-*w(rn(fc)+rf(fc))  e-it(vn(k)+vj(k))  ^  ^ 


is  the  (phase  compensated)  cross  ambiguity  function  be¬ 
tween  the  received  signal  on  antenna  n  and  the  transmit¬ 
ted  signal  on  antenna  j  over  block  k.  Here  the  delay  and 
doppler  terms  Tn(k)  and  vn(k)  are  related  to  the  state  vari¬ 
ables  9  via  (3).  Substituting  in  (6)  yields 
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2NEk 


N- 1 


E 


n,j=0 


(9) 


The  approximate  MLEs  for  the  delay-doppler  and  reflection 
parameters  2  are  obtained  by  maximising  this  function,  ie 


9(k) 

a(k ) 


argmax060(fc) 


N-l 

E  *!>> 


n,j= 0 


1 

~NE~k 


N-l 

E  x(n]m)) 


n,j= 0 


(10) 


xThe  cross  terms  become  insignificant  as  N  — »  oo. 

2 We  don’t  use  the  reflectivity  phase  information  here. 


3 


SNR  (dB) 

Bias 

Std.  Dev. 

r 

r 

P 

P 

r 

r 

P 

P 

oo 

0.14 

41 

0 

0 

0 

0 

0 

0 

40 

-15 

0.02 

0 

0 

74 

6.2 

0 

0 

30 

-28 

31 

0 

0 

175 

14 

0 

0 

20 

-92 

-0.15 

0 

0 

350 

3.9 

0 

0 

10 

-193 

6.6 

0 

0 

381 

16 

0 

0 

TABLE  I 

Bias  and  Std.  Dev.  of  Estimates  vs  SNR 


Pos’n 

Bias 

Std. 

Dev. 

r 

r 

P 

P 

r 

r 

P 

P 

1 

-92 

-0.15 

0 

0 

350 

3.9 

0 

0 

2 

-178 

-3.8 

0 

0.1 

119 

2.8 

0 

0.2 

3 

-120 

-0.9 

-1.5 

0 

226 

3.0 

1.1 

0 

TABLE  II 

Bias  and  Variance  of  Estimates  for  3  target  positions 


where  @(fc)  C  R4  is  the  so-called  validation  gate.  The 
statistic  (9)  represents  the  magnitude  of  the  coherent  sum 
of  N2  matched  filters.  As  we  shall  later  see,  neglecting  the 
cross  terms  in  the  likelihood  function  results  in  estimates 
which  can  be  significantly  biased,  especially  for  small  ar¬ 
rays.  A  correction  scheme  for  reducing  the  bias  in  super¬ 
imposed  signal  problems  such  as  the  one  presently  at  hand 
is  described  in  [4]. 

IV.  Simulation  Examples 

A  single  target  starts  at  the  point  in  the  x  —  z  plane  at 
r  =  10  km,  and  </>  =  7t/4,  and  moves  at  speed  167  m/s  to¬ 
wards  the  origin.  At  a  distance  of  1  km  from  the  origin,  the 
target  turns  right  and  heads  in  the  y  direction  at  a  speed  of 
200  m/s.  Our  first  3  examples  consider  the  target  param¬ 
eter  estimation  at  three  points  in  the  trajectory  ;  one  near 
the  start,  one  at  the  point  of  turning,  and  one  near  the  end. 
For  these  examples,  we  chose  an  array  of  3  elements,  spaced 
one-half  wavelength.  The  carrier  frequency  was  5  GHz, 
and  pulse  duration  was  500  /rsec.  Identical  raised-cosine 
pulse  shaping  was  used  for  each  antenna.  The  temporal 
block  length  was  8  pulses,  and  the  pulse  repetition  interval 
(PRI)  was  67  msec.  The  target  reflectivity  had  unit  vari¬ 
ance.  We  show  in  table  I,  the  bias  and  standard  deviation 
estimates  for  the  four  state  parameters  for  several  values 
of  signal-to-noise  ratio  (SNR)  3  for  the  target  at  position  1. 

Table  II  shows  the  bias  and  variance  of  the  estimates  for 
the  3  target  positions  at  SNR  =  20  dB. 

In  order  to  assess  the  effect  of  array  length  on  the  sta¬ 
tistical  properties  of  the  estimates,  we  repeated  the  ex¬ 
periments  for  various  array  lengths  at  an  SNR  of  20  dB. 
In  all  cases,  the  transmitted  code  was  a  random  unitary 

‘!We  define  SNR  as  the  ratio  of  target  reflectivity  variance  to  noise 
variance 
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Bias 

Std.  Dev. 
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P 
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r 

P 

P 

3 

-92 

-0.15 

0 

0 

350 

3.9 

0 

0 

7 

-67 

1.7 

0 

0 

265 

5.8 

0 

0 

13 

-18 

0.2 

0 

0 

143 

4.7 

0 

0 

TABLE  III 

Bias  and  Variance  of  Estimates  for  different  array  sizes 


code,  held  constant  for  each  of  100  independent  realisa¬ 
tions  which  were  averaged  to  obtain  the  results  presented 
in  table  III.  For  array  lengths  of  3  and  7,  the  block  length 
was  8,  whilst  for  array  length  of  13,  the  block  length  was 
14  pulses.  Target  position  1  was  used. 

Figure  3  shows  the  ambiguity  function  for  the  transmit¬ 
ted  signal  on  antenna  1. 

Ambiguity  surface  for  Coded  Pulse 


1.4 

1.2  -v. . 


Fig.  3.  Coded  Signal  Ambiguity  Function 


The  second  experiment  deals  with  estimation  of  the  tar¬ 
get  parameters  over  the  complete  trajectory.  Here,  a  new 
and  randomly  chosen  unitary  code  was  used  for  each  signal 
block.  Figure  4  shows  the  true  and  estimated  target  range 
over  the  trajectory.  Note  that  the  range  estimates  are  the 
raw  estimates  from  the  generalised  MF  ;  there  is  no  post¬ 
detection  tracking  used.  Figure  5  shows  the  radial  velocity 
estimates.  The  SNR  was  20  dB  and  we  used  3  elements 
and  a  block  length  of  8. 

A.  Discussion 

The  above  results  show  that  for  small  array  sizes,  the 
target  parameter  estimates  are  generally  biased.  This  is 
due  to  the  fact  that  certain  cross  terms  are  neglected  in 
the  formulation  of  the  estimator.  These  terms  are  sig¬ 
nificant  for  small  arrays.  The  results  shown  in  table  III 
however  suggest  that  the  bias  problem  may  diminish  as 
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Fig.  4.  Range  Estimation  (3  elements,  SNR  =  20  dB) 


Fig.  5.  Radial  Velocity  Estimation  (3  elements,  SNR  =  20  dB) 


the  array  size  becomes  larger.  All  experiments  show  that 
the  range  estimates  generally  improve  for  higher  SNR  and 
larger  arrays,  but  the  velocity  estimation  is  problematic. 
The  problem  with  estimation  of  the  velocity  (and  some¬ 
times  the  wavenumber  parameters)  is  due  to  a  combination 
of  bias  and  the  presence  of  local  maxima  in  the  estimation 
statistic.  This  behaviour  is  particularly  evident  in  figure  5, 
where  the  lack  of  effective  doppler  gating  results  in  the  pa¬ 
rameter  estimates  moving  to  an  incorrect  local  maximum. 
All  these  results  should  be  considered  with  some  degree  of 
caution  as  we  have  made  no  attempt  to  optimise  the  trans¬ 
mitted  codes,  in  particular,  as  such,  no  beamforming  has 
been  performed. 

V.  Conclusions  and  Further  Work 

In  this  paper,  we  have  proposed  a  new  approximate  max¬ 
imum  likelihood  estimator  for  the  parameters  of  a  single 
target  when  a  multiple  antenna  radar  is  used.  The  pro¬ 


posed  estimator  combines  both  the  usual  delay-doppler  es¬ 
timation  with  array  beamforming,  by  the  use  of  multiple 
coded  transmitter  waveforms.  We  have  examined  the  per¬ 
formance  of  the  estimator  using  random  unitary  transmit¬ 
ter  codes.  The  performance  of  the  receiver  in  the  absence  of 
any  attempt  to  optimise  the  codes  is  problematic,  however 
results  to  be  presented  in  [6],  show  significant  improvement, 
mainly  due  to  the  inherent  use  of  near-optimal  beamform¬ 
ing  as  a  result  of  the  optimisation  of  transmitter  code  se¬ 
lection. 

In  future  work  [6],  we  derive  an  approximate  Cramer- 
Rao  bound  for  the  receiver  presented  here,  and  use  this 
bound  to  optimise  the  choice  of  transmitter  codes.  We  also 
investigate  incorporating  bias  correction  using  the  scheme 
proposed  in  [4].  The  paper  [7]  is  concerned  with  optimising 
tracker  performance  when  a  Kalman  filter  based  tracker 
is  used  to  process  the  raw  parameter  estimated  from  the 
approximate  ML  receiver. 
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Abstract —  This  paper  addresses  the  problem  of  waveform 
design  for  general  diversity  systems.  The  paper  firstly  in¬ 
troduces  a  general  model  for  such  systems,  and  then  con¬ 
siders  linear  and  nonlinear  cases.  Examples  of  each  case 
are  given  -  in  the  linear  case,  the  MIMO  communications 
design  problem  ;  in  the  nonlinear  case,  the  MIMO  radar 
waveform  design  problem.  In  the  latter  case,  a  simulation 
example  is  provided  which  illustrates  the  potential  benefits 
which  might  be  obtained.  Finally,  the  paper  concludes  with 
a  brief  discussion  of  the  MIMO  radar  tracking  problem. 

I.  Introduction 

There  has  been  considerable  interest  in  the  last  several 
years  among  wireless  telecommunications  researchers  in  the 
use  of  diversity  techniques  to  improve  the  capacity  and  per¬ 
formance  of  such  systems.  In  particular,  the  use  of  multi¬ 
ple  transmit  and  receive  antennae  has  been  demonstrated 
to  offer  significant  performance  improvement  in  wireless 
systems.  More  recently,  the  idea  of  exploiting  diversity  in 
radar  systems  has  also  drawn  attention.  At  this  conference, 
there  were  a  number  of  important  papers  presented  which 
focused  on  various  fundamental  aspects  of  Multiple-Input, 
Multiple-Output  (MIMO)  radar  systems.  Also,  a  confer¬ 
ence  dedicated  to  radar  diversity  and  waveform  design  was 
recently  held  in  Edinburgh,  UK,  where  many  aspects  of 
MIMO  radar  were  presented.  We  can  expect  that  the  area 
of  MIMO  radar  will  receive  an  increasing  amount  of  atten¬ 
tion  by  researchers  and  research  sponsors. 

The  purpose  of  the  recent  work  in  the  area  of  MIMO 
radar  being  conducted  at  the  University  of  Adelaide  is  four¬ 
fold  :  (i)  Firstly  we  seek  to  understand  how  models  used 
in  MIMO  communications  systems  and  associated  trans¬ 
mitter  and  receiver  design  methods,  might  be  relevant  to 
the  radar  case,  (ii)  Secondly,  we  seek  to  study  how  to 
design  appropriate  receiver  structures  for  both  the  radar 
and  communications  cases,  and  importantly,  be  able  to  of¬ 
fer  some  guidance  as  to  the  statistical  performance  of  such 
receivers,  (iii)  Thirdly,  we  seek  to  address  the  transmit¬ 
ter  waveform  design  problem  in  each  case.  We  envisage  a 
system  whereby  the  transmitted  waveforms  are  adaptively 
derived  from  current  system  state  information  under  some 
performance  criterion  eg  tracking  MSE.  (iv)  Fourthly,  we 
are  interested  in  how  feedback  information  is  to  be  sent 
back  to  the  transmitter  in  the  case  when  such  information 
must  be  sent  over  capacity  limited  channels. 

In  [1],  we  considered  channel  modelling  and  receiver 
design  for  MIMO  communications.  The  proposed  re¬ 


ceiver  used  one-step  predictions  for  the  channel  gains  in¬ 
corporated  into  a  Per-Survivor  Processing  based  detection 
scheme.  Both  [2]  and  [3]  dealt  with  MIMO  transmit  code 
design  for  the  communications  problem.  In  these  papers, 
communications  based  criteria  such  as  maximal  distance 
between  codes  and  maximal  diversity  were  applied.  In  [4], 
channel  tracking  MSE  was  used  as  a  code  design  criterion, 
and  was  observed  to  offer  better  performance  at  low  SNRs. 
It  was  demonstrated  in  [5]  that  channel  estimation  MSE  is 
not  the  optimal  criterion  from  the  point  of  view  of  sym¬ 
bol  detection.  We  started  to  address  the  MIMO  radar 
problem  in  [6]  where  channel  modelling  and  receiver  de¬ 
sign  were  considered  in  the  single  target  case.  Importantly, 
the  receiver  structure  considered  incorporated  a  number  of 
constraints  on  the  transmit  codes  in  order  to  permit  the 
calculation  of  the  Cramer-Rao  bound  on  target  parameter 
estimation  performance.  This  is  required  for  parameterisa- 
tion  of  the  tracker  following  the  parameter  estimation  step. 

The  purpose  of  this  paper  is  to  provide  a  general  descrip¬ 
tion  of  the  MIMO  code  design  problem,  and  to  consider  two 
special  cases.  The  first  is  the  case  where  the  model  is  linear 
in  the  target /channel  parameters  and  the  MSE  estimation 
performance  is  time-invariant.  An  example  is  the  commu¬ 
nications  problem  as  discussed  in  [4] .  The  second  example 
is  the  MIMO  radar  problem  as  considered  in  [6]  in  which 
the  target  parameters  appear  in  the  model  in  a  nonlinear 
way.  Here  we  extend  that  work  to  specify  how  the  codes 
may  be  adaptively  chosen  to  optimise  the  target  parame¬ 
ter  estimation.  We  then  mention  how  the  tracking  problem 
is  set  up,  with  results  to  be  provided  in  work  to  be  pub¬ 
lished  [7].  Finally,  we  conclude  with  some  open  research 
questions. 

II.  MIMO  System  Models 

We  assume  a  channel/target  model  with  linear  dynamics, 
and  generally,  non-linear  measurements  of  the  form 

xk+i  =  A  xk  +  B  uk 

Vk  =  fk(xk)  +  vk,  (1) 

where  xk  £  CN  is  the  target /channel  state,  uk  £  Cp  is 
the  white  state  noise  process,  yk  £  CM  is  the  white  mea¬ 
surement  noise  process,  and  vk  £  CM  denotes  the  measure¬ 
ment  noise  process.  We  assume  that  the  matrix  A  is  stable, 
and  that  the  noise  processes  uk  and  vk  are  zero  mean  Gaus¬ 
sian  with  known  covariance  matrices  Q  and  R  respectively. 


The  matrices  A  and  B  specify  the  target /channel  dynamics 
whilst  the  functions  fk(-)  describes  the  mapping  from  the 
target  state  to  the  observed  measurements. 

A.  The  Linear  Case 

In  the  linear  case,  fk{ x)  =  Ck  Xk,  where  Ck  are  M  x  N 
matrices.  An  example  of  the  linear  case  is  in  the  MIMO 
communications  problem,  [4]  with  n  transmit  antennae  and 
£  receive  antennae.  The  states  Xk  are  the  complex  gains  1 
associated  with  a  Rayleigh/Ricean  fading  channel  with  n 
inputs  and  t  outputs,  and  have  dimension  N  =  n£.  The 
measurement  vector  yk  has  length  M  =  mi  where  m  is 
the  temporal  block  length.  The  matrix  Ck  has  the  form 
Ck  =  It®  Xk  P  where  Xk  is  a  sequence  of  unitary  m  x  n 
matrices  which  are  chosen  from  a  set  U  of  cardinality  2P  in 
order  to  encode  p  input  message  bits  per  temporal  block. 
Here  P  is  a  fixed  precoder  matrix  of  size  nxn.  Our  design 
task  is  to  choose  M  based  on  the  a  priori  statistics  of  the 
channel,  in  order  to  minimise  in  an  appropriate  sense,  the 
channel  estimation  error.  The  MLE  of  the  state  term  Ck  Xk 
given  yk  is  simply  yk  itself,  so  in  reference  to  the  block 
diagram  (figure  1),  the  MLE  block  is  trivial.  We  return  to 
this  problem  in  section  III. 


III.  MIMO  Communications  Example 

In  a  communications  problem,  the  overall  code  design 
objective  is  to  minimise  the  resulting  symbol  detection  er¬ 
ror  probability.  This  is  an  unsolved  problem  in  general, 
although  results  exist  for  the  temporally  uncorrelated  and 
asymptotically  high  SNR  case  (eg  [9],  [10]).  Here  we  take 
the  approach  of  minimising  the  channel  estimation  mean- 
square  error,  although  it  is  known  that  this  does  not  gen¬ 
erally  lead  to  optimal  symbol  detection  performance  [5]. 
However  the  codes  so  designed  can  offer  improvement  over 
the  existing  designs  in  some  cases,  eg  at  low  SNR.  The 
optimal  tracker  is  given  by  the  Kalman  predictor  (causal¬ 
ity  is  generally  required  for  any  receiver  incorporating  the 
channel  estimates) 

Xk+l\k  —  A  ( I  GkCk)  Xk\k  —  1  T  Gk  yk  (3) 

where  the  gain  Gk  is  given  by 

Gk  =  EfcCf  {Ck^kCl  +  R)~ 1  ,  (4) 

where  E k  satisfies  the  Riccati  equation 


B.  The  Nonlinear  Case 

In  the  nonlinear  case,  we  construct  an  approximate  MLE 
of  Xk  given  each  measurement  yk-  Such  an  estimate  ne¬ 
glects  the  a  priori  dynamics  of  the  state,  these  being  taken 
into  account  in  the  tracker.  We  will  assume  that  these  esti¬ 
mates  are  given  by  Xk  =  g(yk)  for  some  function  g(.)  which 
is  generally  nonlinear.  For  example,  g(.)  might  represent 
a  generalised  matched-filter  correlator  such  as  will  be  de¬ 
scribed  in  the  MIMO  radar  example  in  section  IV.  Follow¬ 
ing  the  idea  of  Kershaw  and  Evans  [8],  we  characterise  the 
statistical  performance  of  the  estimator  by  assuming  that 
it  is  unbiased  and  efficient,  so  that  Xk  =  Xk  +  Wk  where 
E{wfc}  =  0  and  E{u> kXu^}  =  F_1  where  F  is  the  Fisher 
information  given  by 


F=  - 


d2p(yk\xk) 

dxkdxl 


(2) 


This  simplifying  assumption  is  made  because  it  is  gen¬ 
erally  not  possible  to  determine  the  exact  statistical  per¬ 
formance  for  any  exact  or  approximate  MLE  in  the  non¬ 
linear  case.  In  general,  F  depends  on  the  state  Xk  as  well 
as  the  underlying  code  Ck  (specified  in  the  form  of  /*,(.)). 
One  possible  objective  might  be  to  choose  the  code  to 
maximise  the  objective  function  J(Ck,Xk )  =  Tr (FWF)  = 
\\FW1/2\\2f,  where  IF  is  a  diagonal,  state-dependent  ma¬ 
trix.  Of  course,  since  we  don’t  know  Xk  we  need  to  replace 
its  value  by  an  estimate  (perhaps  from  the  tracker)  and  per¬ 
form  the  optimisation  in  a  block  by  block  adaptive  manner. 
We  return  to  this  problem  in  section  IV 


[  n  the  following,  we  use  real  variables  for  simplicity  reasons,  so  the 
model  for  the  complex  case  is  built  from  the  real  and  imaginary  parts 


Sfc+1  =  A  (r,k  -  EfcCj  (OfeEfcCj  +  R)  1Ck^k)AT 
+BQBt  .  (5) 

We  can  show  that  Tr(Efc)  is  asymptotically  constant,  and 
is  independent  of  the  choice  of  code  Xk  and  only  depends 
on  the  (fixed)  precoder  P.  Indeed,  Tr(Efc)  — »  Tr(E)  where 
E  is  the  solution  to  the  algebraic  Riccati  equation 


E  =  A  (S  -  E PT  (PSPT  +  R)  1  Ps)  At  +  BQBt  .  (6) 

In  [4],  we  show  how  to  compute  the  Jacobian  (first 
derivative)  of  E  with  respect  to  P.  We  then  define  a  gradi¬ 
ent  projection  scheme  for  minimising  Tr(E)  over  P  subject 
to  the  power  constraint  Tr(PPT)  =  1.  These  results  rely 
on  methods  of  matrix  differential  calculus  and  perturbation 
theorems  for  eigenvectors.  The  reader  is  referred  to  [4]  for 
details  and  some  simulation  results. 

IV.  MIMO  Radar  Example 

In  this  section,  we  describe  the  signal  model  and  receiver 
for  a  MIMO  single  target  problem.  In  the  MIMO  radar 
problem,  the  baseband  signal  transmitted  on  antenna  j  for 
temporal  block  k  has  the  form 


m—  1 

Sj,k{t)  =  XiAk)  9j(t  -  (km  +  q)T)  ,  (7) 

9=0 

where  the  gj(-)  are  unit  energy  pulse  shaping  functions 
supported  on  [0,  T],  and  X(k)  is  the  space-time  code  sent 
for  block  k.  The  baseband  received  signal  on  antenna  r  is 


n—  1 

yr,k(t)  =  a(k)J2eitiMk)+Vim  SjAt-TrW-Tjik)) 

3=0 

xe^-i  +wrJ(k)  .  (8) 

Here  Tj(k)  and  Vj(k)  are  the  one-way  delay  and  doppler 
between  antenna  j  and  the  target  assumed  constant  over 
the  duration  of  each  block,  and  a(k)  is  the  magnitude  of  the 
target  reflectivity  coefficient  together  with  the  path  loss. 
The  terms  4>rj  are  uniform  iid  variables  on  [ — 7r,  7r] .  The 
delay  and  doppler  parameters  contained  in  (8)  are  functions 
of  the  delay  and  doppler  respectively  to  some  reference  an¬ 
tenna,  together  with  the  array  geometry  and  target  motion. 
More  precisely,  for  a  uniform  linear  array  with  spacing  d 
metres,  to  linear  terms,  2 

Tj  =  t0  —  jdp/c 

Vj  =  v0-  jdcop/c  ,  (9) 

where  To  =  r/c  and  vq  =  fco/c  are  the  reference  delay  and 
doppler,  p  is  the  spatial  wavenumber,  and  p  its  time  deriva¬ 
tive.  Here  r  and  r  are  the  radial  distance  from  the  reference 
antenna  to  the  target  and  its  time  derivative  respectively, 
io  is  the  carrier  frequency  (rad/s),  and  c  is  the  propagation 
speed.  We  will  denote  the  state  of  the  system  (target)  by 
the  vector  xk  =  [a(k),  r(k),  r(k),  p(k),  p{k)]T .  Notice 
that  (8)  is  a  nonlinear  mapping  from  the  target  state  to 
the  observed  signal. 

A.  Approximate  MLE 

In  [6],  we  described  a  generalised  matched  filter  receiver 
for  this  MIMO  radar  model.  The  data  likelihood  con¬ 
tains  cross-terms  between  the  individual  receive  antenna 
signal  models  which  make  exact  MLEs  impractical  to  ob¬ 
tain.  These  cross  terms  are  known  to  tend  to  zero  as  the 
data  length  goes  to  infinity.  However,  the  receiver  proposed 
in  [6]  imposes  a  number  of  constraints  on  the  transmitter 
MIMO  codes,  among  which  is  that  they  be  unitary.  This  re¬ 
sults  in  the  suppression  of  the  aforementioned  cross  terms, 
and  thus  this  receiver  is  a  better  approximation  to  the  ML 
one. 

More  specifically,  the  data  log  likelihood  for  a  block  of 
received  data  is  given  by 

JV-l  N—l 

t  ~a  X!  si  -  T"  “  Ti)  e*071'" 

**  n— 0  j= 0 

x  £**("» +".#)  2  dt  .  (10) 

Neglecting  the  cross  terms,  we  use  instead  the  approxi¬ 
mate  log  likelihood 

2In  the  following,  we  drop  the  time  block  index  k 


N2a2 

l=a  ^2  Re  {e-*0"’JXn,fc(T„  +  rk,  vn  +  vk)} - — 

n,j= 0 

(11) 

which  is  maximised  by  maximising  the  detection  statistic 

N-l 

V(r,  r,  p,  p)  =  I Xj,k  [(t  -  (. j  +  k)dp)/c, 

j,k=0 

(r  -  (j  +  k)dp)u/c]\2  ,  (12) 

where 

Xi,fc(r,  v)  =  J  rj(t )  s%(t  -  r)  e~ltv  dt  ,  (13) 

are  the  outputs  of  the  generalised  matched  filters  for  each 
receive  antenna  rj,  corresponding  to  each  transmitted  sig¬ 
nal  sk.  Approximate  ML  estimates  for  r,  f,  p  and  p  are 
obtained  by  maximising  (12),  and  an  estimate  of  a  is  ob¬ 
tained  from  that  maximum  value  of  rj.  Phase  estimates  can 
also  be  obtained  but  are  not  used  here. 

In  [6],  we  were  able  to  evaluate  the  Fisher  information 
matrix  corresponding  to  (8)  by  imposing  a  number  of  con¬ 
straints  on  the  signal  codes.  Indeed,  these  constraints  have 
the  form  of  moments 

XXH  =  I 

XT  jXH  =  A  j  ,  (14) 

for  j  =  1,  2,  where  the  Tj  are  specified  diagonal,  posi¬ 
tive  matrices  which  depend  on  the  pulse  shaping  functions, 
and  the  A  j  are  free  diagonal  (positive)  matrices.  In  this 
case,  the  Fisher  information  F  is  independent  of  the  target 
state,  and  only  depend  on  the  choice  of  code  X  via  the 
free  diagonal  terms  of  Aj  and  A2.  Thus  the  optimisation 
task  is  to  maximise  J(Ai,  A^,x)  =  ||F(Ai,  A2)  W1^2{x)\\2F, 
where  IF  is  a  possibly  state  (ie  x)  dependent  positive  di¬ 
agonal  weighting  subject  to  the  constraints  (14).  We  used 
a  projection  gradient  approach,  but  the  issues  concerning 
the  feasibility  of  the  constraints  (14)  and  the  convergence 
of  the  scheme  remain  open  questions  which  we  aim  to  ad¬ 
dress  in  [7]. 

We  simulated  the  system  described,  with  the  target  lo¬ 
cated  10km  from  the  antenna  array,  moving  towards  the 
array  at  speed  267  m/s,  at  wavenumber  0.1,  and  wavenum¬ 
ber  derivative  0.01.  The  results  shown  in  table  1  show  the 
proportional  errors  obtained  over  50  simulation  runs  each 
consisting  of  100  pulse  blocks.  We  used  3  antennas  and  a 
temporal  block  length  of  4  pulses.  Raised  cosine  shaping 
with  pulse  duration  50  msec  was  used.  The  SNR  was  20  dB 
(total  energy  per  pulse  divided  by  noise  variance  per  sen¬ 
sor).  We  compared  2  code  sets  -  one  random  unitary  codes, 


and  the  other  optimal  codes  chosen  according  to  the  above 
criterion  with  weighting  matrix  W  =  diag(0.01, 1, 10, 10). 
A  significant  improvement  is  noted  when  the  optimal  codes 
are  used. 


Parameter 

range 

radial  vel. 

P 

P 

Opt  bias  (%) 

2.5 

0.02 

0.6 

1.5 

Opt  stdev  (%) 

1.7 

0.004 

1.4 

1.5 

Ran  bias  (%) 

4.3 

0.07 

0.8 

2.5 

Ran  stdev  (%) 

3.7 

0.03 

2.6 

2.2 

TABLE  I 

Bias  and  Standard  Deviation  Estimates  for  Optimal  and 
Random  Codes 


B.  Discussion 

As  shown  by  Furhmann  [11],  the  use  of  unitary  codes 
results  in  no  directivity  in  the  transmitted  energy  which 
whilst  desirable  initially  in  the  absence  of  a  priori  infor¬ 
mation  about  the  target  state,  is  clearly  sub-optimal  when 
estimates  of  the  target  state  are  available.  However,  we 
found  it  necessary  to  impose  a  number  of  simplifying  con¬ 
straints  on  the  codes,  including  the  unitary  one,  in  order  to 
evaluate,  even  approximately,  the  statistical  performance 
of  the  (approximate)  MLE.  More  work  is  thus  required  in 
the  area  of  receiver  performance  evaluation  so  that  the  uni¬ 
tary  constraints  might  be  relaxed.  Alternatively,  we  could 
explicitly  add  beamforming  into  the  code  design  problem. 
These  are  areas  for  ongoing  research. 

V.  Signal  Design  for  Tracking 

For  tracking  problems,  we  feed  the  “raw”  estimates  of 
the  target  state  xk  to  a  Kalman  filter  derived  from  the 
state  space  model 


Xk+i  =  A  xk  +  B  uk 
Xk  =  xk  +  wk 


(15) 


where  wk  ~  A/"(0,  F1-1).  The  resulting  tracking  error  co- 
variance  is  given  (approximately)  by  the  Riccati  difference 
equation 


^fc+i|fc+i 


AZk\kAT  +  BQBt 

^Wi|fc  —  5]fc+i|fe  (Sfc+i|fe  +  F  X)  5]fc+i|fc 

(16) 


where  F  depends  on  the  chosen  code  as  described  in  the 
previous  section.  Thus  we  seek  to  select  Xk+\  by  minimis¬ 
ing  Tr{Efc+1|k+1}  with  xk+i  replaced  by  the  current  one 
step  prediction  state  estimate  xk+i\k-  Matrix  differential 
calculus  can  be  used  to  obtain  the  Jacobian  of  Tr(Efc+1|fc+1) 
wrt  the  Ai  and  A2  matrices,  and  an  alternating  convex  pro¬ 
jection  gradient  algorithm  thus  derived.  We  describe  this 
approach  in  [7]. 


VI.  Conclusions  and  Future  Research 

In  this  paper,  we  have  introduced  a  general  model  for 
waveform  design  in  MIMO  diversity  systems.  We  have 
considered  both  linear  and  non-linear  models,  and  have 
indicated  methodologies  for  optimal  code  design  in  each 
case.  Examples  of  designs  for  MIMO  communications  sys¬ 
tems  (linear),  and  for  MIMO  radar  systems  (nonlinear) 
have  been  given.  In  the  latter  case,  we  have  included  a 
simulation  example  illustrating  the  potential  benefit  of  de¬ 
signing  the  transmitter  code  to  optimise  the  performance 
of  an  approximate  ML  receiver.  Finally  we  have  indicated 
how  to  design  optimal  codes  for  tracking  systems. 

Several  open  questions  remain,  including  receiver  perfor¬ 
mance  characterisation  for  non-unitary  codes,  convergence 
of  numerical  schemes  for  code  optimisation,  using  codes 
chosen  from  discrete  sets,  performance  analysis  of  tracking 
systems  and  the  implementation  of  feedback  information 
in  the  system. 
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Abstract — This  paper  addresses  the  problem  of  transmitter 
design  for  a  single  target  MIMO  radar  system.  The  transmitted 
space-time  codes  are  unitary  apart  from  the  inclusion  of  a 
transmit  beamformer.  The  paper  demonstrates  that  the  use 
of  unitary  codes  leads  to  a  simple  coherent  correlations-based 
maximum-likelihood  receiver  yielding  unbiased  estimates  of  the 
target  range  and  wavenumber,  and  their  time  derivatives.  The 
transmitter  codes  are  chosen  to  maximise  the  Fisher  information 
for  the  model.  This  leads  to  an  alternating  projection  gradient 
adaptive  algorithm  for  dynamic  code  design.  The  paper  argues 
that  proper  choice  of  beamformer  phases  can  lead  to  efficient 
estimators,  thus  validating  the  optimisation  methodology. 

I.  Introduction 

In  recent  years,  there  has  been  a  rapid  increase  in  interest  in 
the  use  of  multiple  input  -  multiple  output  (MIMO)  systems 
for  radar  detection,  estimation  and  tracking.  Our  recent 
work  has  focused  on  receiver  design  [1],  and  transmitter 
design  [2].  In  each  of  these  papers,  the  transmitter  codes 
were  constrained  to  be  orthogonal,  leading  to  a  simple 
derivation  of  the  maximum  likelihood  (ML)  receiver. 
Additional  orthogonality  constraints  were  included  to  permit 
the  derivation  of  the  corresponding  Fisher  information  for  the 
model.  Following  the  methodology  proposed  by  Kershaw  and 
Evans  [7]  in  the  SISO  case,  we  attempt  to  optimise  system 
performance  by  maximising  the  Fisher  information  of  the 
model  by  appropriate  choice  of  codes.  In  the  SISO  case, 
the  ML  estimators  for  range  and  Doppler  are  efficient  for 
finite  sample  sizes,  thus  validating  this  approach.  However,  it 
is  unclear  whether  this  property  holds  for  the  MIMO  case, 
where  the  time  derivatives  of  range  and  Doppler  are  also 
estimated.  Another  important  issue  is  that  the  use  of  unitary 
codes  precludes  directivity  being  introduced  in  the  transmitter 
beamformer  [6],  This  is  obviously  a  clear  deficiency  in  a 
MIMO  design  approach. 

In  this  paper,  we  include  explicit  (linear)  transmit 
beamforming  in  addition  to  unitary  space-time  coding. 
However,  retaining  orthogonality  between  transmitted 
waveforms  (so  that  a  simple  ML  receiver  can  be  designed), 
means  the  resulting  Fisher  information  is  independent  of  the 
beamformer  phases.  Thus  it  is  unlikely,  in  general,  that  the 
resulting  estimates  are  efficient.  A  natural  question  arises 
as  to  whether  efficiency  can  be  attained  by  “steering”  the 
transmitter  beam  towards  the  expected  target  location.  In  a 


more  general  tracking  context  [3],  the  tracker  can  assist  with 
this  operation. 

The  layout  of  the  paper  is  as  follows.  Firstly  we  introduce 
the  coherent  signal  model  for  single  target  MIMO  radar,  and 
thus  derive  a  correlation  based  ML  receiver.  We  then  provide 
a  characterisation  of  the  receiver  performance  in  terms  of  the 
transmit  codes  (including  beamformer).  We  then  introduce  an 
procedure  for  the  design  of  optimal  codes.  Simulation  results 
will  be  presented  at  the  meeting,  and  are  also  available  at  [4], 

II.  Signal  Model 

We  consider  a  radar  system  with  N  transmit  and  N  receive 
antennae.  The  baseband  transmitted  signal  for  temporal  block 
k  and  transmit  antenna  n  is  given  by 

M—l 

sn,k{t)  =  ^2  Yn,m(k)  gn(t  -  ( kM  +  m)T)  ,  (1) 

m= 0 

where  M  >  N  is  the  number  of  pulses  in  each  temporal 
block,  T  is  the  pulse  duration  and  gn  is  the  baseband  pulse 
shaping  function  for  antenna  n,  which  is  assumed  real  and 
twice  differentiable,  supported  on  [0,  T],  and  with  gn( 0)  = 
gn(T)  =  0.  The  quantity  M  denotes  the  pulse  repetition 
interval  (PRI).  Here,  Y{k)  £  <CNxM  is  the  transmitted  space- 
time  code  for  block  fc,  given  by 

Yn,m(k)  =  Wn{k)  Xnm{k)  ,  (2) 

where  X(k)  £  £,NxM  is  an  unitary  code  (ie  XXH  =  7y), 
and  W„(k)  denotes  the  transmitter  beamformer  weights 
for  block  k.  The  signals  are  modulated  onto  a  carrier  with 
frequency  /o  (Hz). 

We  assume  a  single  target  located  at  a  radial  distance  r  with 
respect  to  the  origin  located  at  a  reference  element  of  the 
transmit  array.  For  simplicity  of  presentation,  we  assume  in  the 
following,  that  the  transmit  and  receive  arrays  are  uniformly 
spaced  linear  arrays  located  along  the  y  axis,  and  are  co¬ 
located.  The  case  where  there  is  a  constant  distance  between 
the  arrays  is  easily  dealt  with.  We  address  the  case  where 
transmit  and  receive  arrays  are  in  relative  motion  in  [5].  To 


linear  terms,  the  two-way  delay  and  Doppler  from  element  n 
to  element  m  are 


chosen  subset  of  the  4  dimensional  parameter  space.  The  mean 
value  of  the  statistic  is 


7"n,m 

L'n.iti 


(k) 

(k) 


2  r(k)  —  { n  +  m)dp(k) 

(2 r(k)  -  (n  +  m)dp(k))  f0 
c 


(3) 


where  d  is  the  array  spacing,  c  the  propagation  speed,  p 
is  the  target  wavenumber,  with  p  being  its  time  derivative, 
and  r  is  the  radial  target  velocity.  For  the  purposes  of  target 
tracking,  we  shall  define  the  target  state  to  be  the  vector 
6  =  (a,  <p.  r,  r,  p,  p)  which  is  a  function  of  the  temporal  block 
index  k.  The  variables  (a,  tj>)  are  defined  below  and  (r,  r,  p,  p) 
are  defined  in  [1], 


The  basebanded  received  signal  at  antenna  m  for  transmit 
block  k  for  a  fully  coherent  system  1  is 


JV-l 

rm,k(t)  =  a(k)  el<t>(k)  ^  sn.k  ( t  -  Tn>m(k)) 

n=  0 

x  e2 nitvn,m(k)  e-2Kif0Tn,m(k)  i(4) 

where  a(k)  >  0  is  the  target  reflectivity  amplitude,  and  tj>{k) 
denotes  the  associated  phase  shift.  The  quantities  £m4(f) 
represent  receiver  noise,  and  are  assumed  to  be  zero  mean, 
white,  Gaussian  processes  with  identical  variance  o'2,  and  in¬ 
dependent  between  receivers.  The  approximate  log  likelihood 
(neglecting  cross  terms)  is  (dropping  the  k  index  for  clarity) 
[8] 


{JV  — 1 

Xn1m  j 

n,m— (I 


a2N 

2 


,(5) 


where  we  assume  that  [|s„jt||2  =  pnki  with  J2npn,k  =  1  for 
all  /,',  and 


Xn,m.{r,u)  =e2w,foT  j  rm(t)  s*n(t,  -  r)  e  2lvlvi  dt  ,  (6) 

denotes  the  matched  filter  for  transmit  waveform  n  received  on 
antenna  m.  Maximising  (5)  over  a  and  <f>  (assumed  determinis¬ 
tic  and  unknown)  we  obtain  the  sufficient  statistic  (generalised 
matched  filter),  written  as  a  function  of  the  target  state 


V  (r,r,p,p) 


1 

2  No2 


N—l 

^2  Xn,m{T{6),v{6)) 

n,m= 0 


,  (7) 


which  is  an  coherent  sum  of  the  outputs  of  the  N2  matched 
filters.  In  order  to  estimate  the  target  state,  we  perform  a  con¬ 
strained  maximisation  of  (7)  subject  to  the  constraints  given 
in  (3).  This  generally  requires  a  search  over  an  appropriately 


'We  assume  all  transmitters  and  receivers  are  phase  locked  to  a  reference 
clock. 


E  r}(r,  r,  p,  p)  =  +  1/2  .  (8) 

We  thus  observe  a  linear  increase  in  the  mean  detection 
statistic  with  the  number  of  antenna  provided  the  target  signal 
to  noise  ratio  is  sufficiently  large. 

III.  Characterising  Receiver  Performance 

In  earlier  work,  [1]  we  derived  the  Fisher  information  matrix 
for  the  receiver  described  above  when  the  transmit  codes  were 
unitary.  Part  of  this  analysis  also  showed  that  the  cross  terms 
neglected  in  the  derivation  of  (7)  are  suppressed  by  the  use 
of  unitary  codes.  Thus  our  receiver  achieves  close  to  ML 
performance  in  this  case.  It  is  easy  to  demonstrate  that  the 
inclusion  of  the  transmit  beamformer  W  does  not  alter  this 
property.  The  beamformer  is  required  because  no  directivity 
can  be  achieved  (apart  from  that  inherent  in  the  array  elements) 
with  unitary  codes  [6],  As  shown  in  [7],  for  the  SISO  case, 
the  Fisher  information  depends  essentially  on  the  first  two 
time  and  frequency  moments  of  the  transmitted  signal.  In 
order  to  find  tractable  form  for  the  Fisher  information  in  the 
MIMO  case,  we  need  to  make  a  number  of  orthogonality 
assumptions  on  the  transmit  signals  and  their  moments  as  in 
[1].  It  is  a  natural  question  to  ask  as  to  whether  it  is  useful 
to  minimise  a  lower  bound  on  receiver  performance  (Cramer- 
Rao  bound).  Clearly,  if  the  receiver  is  consistent,  then  it  is  a 
sound  approach.  In  the  SISO  case,  it  is  known  [8],  that  the 
matched  filter  (ambiguity  function  peak)  yields  unbiased  and 
consistent  estimators  for  delay  and  Doppler  (if  side  lobes  are 
negligible).  However,  it  is  not  clear  that  maximisation  of  (7) 
provides  consistent  estimators  of  6  for  finite  sample  sizes.  We 
have,  however  demonstrated  that  these  estimators  are  unbiased. 

IV.  Optimising  Receiver  Performance 

We  shall  use  as  our  measure  of  receiver  performance,  the  cost 
function 


■J(Y-8 )  =  log  det  F{Y;8)  ,  (9) 

where  F  is  the  Fisher  information  matrix 


which  is  strictly  positive  definite.  This  quantity  is  a  function 
of  the  quantities  r„  m,  and  thus  of  the  state  variables  r  and 
p,  as  well  as  the  signal  waveforms,  and  thus  of  the  code 
Yn  m.  This  cost  function  is  chosen  because  maximisation  of  J 
implies  minimisation  of  the  Cramer-Rao  lower  bound  on  the 
parameter  estimation  uncertainty.  There  are  four  constraints  to 
be  met.  Firstly,  the  beamformer  weight  magnitudes  specify  the 
transmit  powers  for  each  antenna,  ie  |IE„|2  =  pn  =  g2,  with 
5/  n  q2  =  1.  We  also  have  the  code  constraints 


xk  xn  —  ^n,k 
x'j!  F,  xn  =  0,  V  n  ^  k 
xf  Ttt  x.n  =  0,  V  n  ±  k  ,  (11) 

where  T t  and  1’//  are  diagonal  matrices  specified  by  the 
signalling  pulse,  and  the  x.n  are  the  rows  of  X  (ie  the  sequence 
transmitted  on  antenna  n.  We  can  write 

F(q,  X;  8)  =  ^  q2n  (An{9)  +  B.n{0)m*Ttxn  +  Cn(8)x*Tttx, 

n 

where  An{6),  Bn(6),Cn(6)  are  6  x  6  matrices  whose  entries 
can  be  determined  from  the  form  of  the  second  derivatives  of 
l  wrt  a ,  <j>  and  6.  We  thus  have 


important  assumption  underlying  the  approach  is  that  reduc¬ 
tion  of  the  CRB  leads  to  improved  performance.  Whilst  this 
is  true  in  the  SISO  case,  because  the  delay-Doppler  estimates 
are  efficient  in  that  case,  it  is  not  clear  that  this  property 
holds  in  the  orthogonal  coded  MIMO  case.  We  have  seen  that 
the  CRB  is  independent  of  the  transmit  beamformer  phase, 
suggesting  lack  of  efficiency  of  the  estimates.  Thus  explicit 
transmit  beamforming  is  argued  to  be  essential  in  this  context. 
In  simulation  results  to  be  presented  at  the  meeting  (and 
available  at  [4]),  we  investigate  this  issue  in  detail. 
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Thus  [9] 


d.J(q,  X\6)  =  TrF-1  Vec  dF  ,  (14) 

permits  the  gradient  of  J  wrt  q  and  X  to  be  determined.  We 
can  thus  define  a  projection  gradient  scheme  for  numerical 
constrained  maximisation  of  J  using  the  forms  of  the  projec¬ 
tions  given  in  appendix  III,  applied  to  the  constraints  (13).  Let 
Xj  and  qj  denote  estimates  of  the  parameters  at  iteration  j, 
then  we  perform  the  updates 

Vec  Xj  =  Vec  Xj  +  /j  Tr  /•'  1  ,  Xj ;  6) 

r  dVec  F{qj,Xj-6) 

X  dYec  X 

Qj+i  =  Qj  +  k  Tr  F-1  (qj ,  Xj ;  6) 

„  8Vec  F(qjtXj-,9) 
dq 

Xj+1  =  (V„  Py,  V1J  Xj  ,  (15) 

where  ft  is  a  small  scalar  step  size,  and  9  is  an  estimate  of  the 
target  state  which  is  generally  supplied  by  the  post-receiver 
tracker  [3].  The  projection  steps  are  repeated  f  times  until 
there  is  little  change  in  the  resulting  matrix. 

V.  Conclusion 

This  paper  has  derived  a  simple  maximum  likelihood  receiver 
for  MIMO  radar  systems  with  unitary  space-time  codes  and 
explicit  transmitter  beamforming.  The  transmitter  codes  are 
chosen  to  maximise  the  Fisher  information  associated  with 
the  MIMO  model.  An  gradient  projection  algorithm  has  been 
proposed  to  solve  this  constrained  optimisation  problem.  An 
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Appendix  I 

Derivation  of  Receiver  CRB  -  Outline 


In  order  to  compute  the  Fisher  information  for  the  model  (5), 
we  need  to  evaluate  the  first  and  second  order  derivatives  of 
.Vicml'  ")•  These  are 


dXnjn 

dv 

dXn,rn 

3  t 


d2Xn,r 
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d'2Xn,r 
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d‘2Xr, 


3rd  v 


=  j  rm(t )  °*j  (/  -  r)  dt 

+  2nifoXn,m(T,v) 

=  e27ri/or  J  ^  (/  T)  c  2'""  dt 

/£)«* 

r,n(t)  -gf(t-r)  e-2^  dt 

+2*ifo  ^(r.i/) 

f  rm(f)  t-T )  e~2wivt  dt 
+2nif0dX^m  (t,  v)  .  (16) 


We  then  substitute  for  r,n(t)  from  (4),  and  take  expected 
values.  For  example. 


tv-i 


E  V  (r  V  )  —  C1  e‘4  p27ri/o(r„,m -rfe,m  ) 
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J  Sk(t  -  Tk,m)  S*n(t  -  Tn,m)  dt  .  (17) 
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(20) 


Thus  we  can  determine  the  Fisher  information  terms  : 
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(21) 


We  assume  that  the  differential  time  and  frequency  shifts 
across  the  array  are  negligible  compared  to  the  signal  dura¬ 
tion  and  bandwidth,  and  that  signals  transmitted  on  different 
antennas  are  orthogonal.  Thus 


Tk.m)  *;(/  e~2nit^Vh-m~Vn-m  >  dt 

~  J  |s»(f)|2  dt  Sk,n  , 


(18) 


and  thus 


E  X  ri  .m  {d~n,rm  ^n  .  in  ) 


Pn 


(19) 


where  cg  is  a  constant  depending  only  on  the  pulse  g.  In  (20) 
and  (21), 


v n,m 
(Ski  ^n)t 
^n)^2 


J  tsk(t)  4(f)  dt=(YTtYH)Kn 

j  t2  sk(t)  4(f)  dt=  (YTttYH)k  J22) 


We  can  also  compute  the  mixed  second  order  derivatives  of 
£  which  are  not  shown  here.  We  note  that  all  mixed  second 
order  partials  wrt  a  are  zero. 


We  can  also  show  that  the  expected  value  of  all  first  order 
derivatives  of  £  evaluated  at  the  true  parameter  values,  are 
zero  thus  demonstrating  that  estimates  are  unbiased. 

Appendix  II 

Waveform  Constraints 


This  implies  that  the  expected  value  of  the  test  statistic 

This  process  is  repeated  with  the  terms  in  (16),  assuming  the 
orthogonality  properties  detailed  in  appendix  II  to  give  : 


We  need  to  have  <  sn,sk  >=  Sn  kpn ,  where  pn  >  0  is  the 
energy  transmitted  on  antenna  n  each  block.  This  leads  directly 
to  the  constraint 


(31) 


Y  Yh  =  W  WH  =  A  =  diag  (p0,  pN_k)  ,(23) 

where  we  assume  that  the  baseband  pulses  gn(t)  are  identical 
for  each  n,  are  supported  on  [0,T],  and  have  unit  norm.  Here 
W  is  an  TV  x  N  diagonal  matrix  with  the  array  weights  on  its 
diagonal,  and  A  is  a  diagonal,  positive  definite  matrix  which 
specifies  the  transmit  powers  pn.  Unity  total  transmit  energy 
implies  that  Tr(A)  =  1.  The  second  constraint  is  that 


/f)o* 

t  Sn(t)  dt 

is  zero  whenever  k  /  n.  This  follows  from  the  property 

J  t  g(t  -  IT)  (t-  mT )  dt  =  0  (32) 

for  all  i  ^  m.  When  (.  =  m,  we  have  that 


J  t  sn(t)  s*k{t)  dt  =  qn  5n,k  Y  Tt  Yh  =  At  ,  (24) 

where  L,  is  a  diagonal,  positive  definite  matrix  with  elements 

[Tt]nn  =  J  t\g{t-nT)\ 2  dt  =  j  t\g(t)\2  dt  +  nT  ,  (25) 

and  A,  is  a  diagonal,  positive  definite  matrix.  The  next 
constraint  is 

J  t2  sn(t )  s*k(t)  dt  =  un  SnM  Y  Ttt  Yh  =  Att  ,  (26) 

where  T n  is  a  diagonal,  positive  definite  matrix  with  elements 

P ttU  =  J  t2\g(t-  nT)\2  dt  ,  (27) 

and  An  is  a  diagonal,  positive  definite  matrix. 

The  terms 

jsn(t)d-^dt  =  5>Myn*m 

£,m 

x  J  g(t  —  CT)  ?£(t-mT)  dt  (28) 

are  zero  for  all  k.  n  because  if  g  is  real,  and  g( 0)  =  g(T )  =  0, 
the  integral  on  the  rhs  of  (28)  is  identically  0  for  all  i  and  m. 

Consider  now 


J  t  g(t  -  mT)  ^  (t  -  mT)  dt  =  J  t  g(t)  ^-(t)  dt 

=  -1  ~  J  t  g*(t)  ^(t)  dt  (33) 

leading  to  the  fact  that  the  integral  on  the  lhs  equals  -1/2.  Thus 
this  constraint  is  also  met  by  virtue  of  (23). 

Appendix  III 

Approximation  by  Orthogonal  Matrices 
Suppose  Y  £  (^nxm  wjtj1  f  >  tv  js  given  and  has  full 
rank  N.  We  seek  X  E  CiVxM  satisfying  AT.YH  =  D  where 
T  >  0  is  Hermitian  and  is  given,  and  D  >  0  is  a  free  diagonal 
matrix,  such  that  ||AT  —  Y||  is  minimised.  Let 


Consider  now  the  constraint  XVXH  =  D.  Since  T  >  0,  we 
can  find  V  E  CMxM  such  that  T  =  VVH  with  V  >  0. 
Then  XV  =  A  U  where  A  =  D]  >'2  is  a  free  diagonal  matrix, 
and  U  is  any  N  x  M  matrix  satisfying  UUH  =  I.  Let  U  = 
[  ui  ■  ■  ■  uk  ]  with  each  ut  E  CM.  Thus  X  =  A  U  U_1 
and 

|| Y  -  A||2  =  Tr  (Y  Yh  -  A*  Y  V~H  UH ) 

-  Tr  (A  U  V-1  Yh  +  |A|2  U  r-1  UH)  .  (35) 

We  minimise  this  quantity  over  choice  of  diagonal  A.  We  find 
the  optimal  A  =  diag(<5i  , . . . ,SN)  with 


P  Tff  YH)n,k 

(29)  Then 


[U  V-1  yh]u 

[u  r-1  uH]u  ' 


(36) 


where  Tff  is  a  diagonal  matrix  with  elements 

9  JX,  \\ITV-1  YH]..\2 

||Y  -  A||2  =  Tr  (Y  Yh)  -  ^  1  J"1  ,  (37) 

,  (30)  -'=i  1  ki 

which  we  now  minimise  over  unitary  N  x  M  matrices  U. 
under  the  assumption  that  <y ( 0 )  =  g(T)  =  0.  Notice  that  A  Lagrangian  approach  yields  the  necessary  condition  U  = 
Tff  is  proportional  to  the  identity  matrix,  so  the  required  A  Y  V~H ,  where  A  is  a  free  Hermitian  N  x  .V  matrix, 
orthogonality  constraint  is  met  by  virtue  of  (23).  Thus  I  =  U  U 11  =  A  Y  T_1  Y11  AH ,  or  A-1  A~H  = 

Y  r-1  YH .  We  can  thus  choose  A  =  G-1  where  G  is  the 
unique  symmetric  square  root  of  Y  T_1  YH .  Then  A  = 


-Jsw 


d2g*  dg 

n  -  ( t )  dt.  =  — 
dt 2  v  1  dt 


The  final  constraint  is  that 


A  G~ ]  Y  r_1.  We  can  also  show  that  X  is  unique  and 
is  independent  of  the  choice  of  V .  We  denote  the  mapping 
(projection)  by  X  =  'Pr(5").  Simulations  suggest  that  Y\  is 
a  contraction  mapping  on  KJVxM  with  Frobenius  norm.  In 
the  unitary  case,  where  T  =  D  =  I,  the  optimal  solution  is 
X  =  (Y 1  Y  where  G  is  the  (unique)  symmetric  square  root 
of  YYh. 
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System  Architecture 


Waveform  code  for  next  block  of  pulses 
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Target  Dynamics 
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State  varables:  Q(k)  =  (r(k),r(k),p(k),p(k))T 
with  p{k)  =  sin(9(fc)  sin  (j){k) 

Time  evolution: 


0(fc  +  l)  =  AG(k)  +  Bw(k ) 


where  W(k)  Gaussian  white  noise  process. 


Sub-Pulse 


Coding  of  Pulse 

Technique  of  PulseCoding 
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Transmit  Signal  Model 

Baseband  signal  for  N  array  elements 


M—l 

AQ  =  E  9n(t  ( kM 

m= 0 

where  M  >  N  and  k  pulse-block  number. 

> NxM  . 


m)T ) 


Y(k)  E  C  v  x  M  :  transmitted  space-time  code 


Yn,mi}^)  n(jE)  Xnm{k^ 


Space  -  Time  Coding 

Transmitted  space-time  code  for  pulse-block  k 

Y(k)  e  CNxM 

Indices  -  n  array  element  to  sub-pulse  number  in  fc-th  block 

X(k)  e  CNxM:  Orthogonal  code  ie  XXH  =  ±IN 

Wn  (k):  Array  Beamformer  Weights  for  pulse-block  k 

NB.  Orthogonal  code  enables  directivity  in  Beamformer  in 
contrast  to  Unitary  code. 


-p.  1/1 


Delay  and  Doppler  across  array 


The  two-way  delay  and  Doppler  terms  for  signal 
transmission  from  element  n  to  element  m  after 
reflection  from  a  single  target  are 
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/o  is  the  carrier  frequency 


Received  Signal  at  Array 


The  basebanded  received  signal  at  antenna  m  for 
transmit  block  k  for  a.  fully  coherent  system  fl  is 


N-l 

-  a{k)  el4>{k)  ^  sn,k  (t  -  Tn,m(k)) 

n= 0 

g27TZtZ/n?m  (Al)  ^  27T2/0Tn;m(fc)  _|_  ^ 

a(k)  >  0  target  reflectivity  amplitude  and  <i>(k)  associated  phase  shift. 
£m,k{t)  receiver  noise  Gaussian  processes  with  identical  variance  a2. 

"We  assume  all  transmitters  and  receivers  are  phase  locked  to  a  ref¬ 
erence  clock. 
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Likelihood  Function 

The  method  of  correlation  between  the  transmitted  and  received 
signal  leads  to  the  log  likelihood  function  as 
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where  we  assume 


pn,k  With  Pn,k  =  1  V  k  and 


Xn,m(r,  v)  =  e2mhT  I  rm(t )  s*  (t  -  r)  e  2mut  dt 

which  is  the  generalised  matched  filter. 

This  is  an  approximate  form  since  the  cross  terms  are  neglected. 
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Maximum  Likelihood  Estimate 

Maximisation  of  the  likelihood  function  l  gives  the  sufficient 
statistic  ( the  generalised  matched  flier)  as  a  function  of  the  target 
state  in  the  form 


V(r  ,r,p,p) 
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27V  cr2 


TV— 1 

Y  Xn,m(r(9),iy(9)) 

n,m=0 


This  form  is  a  coherent  sum  of  the  outputs  of  N2  matched  filters. 
We  optimise  rj  subject  to  the  delay  and  Doppler  constraints  for 
the  given  target.  The  mean  value  of  the  statistic  is 

E  rj{r,r,p,p)  =  ^^  +  1/2.  (2) 

Thus  a  linear  increase  of  the  mean  detection  statistic  as  a 
function  of  the  number  of  antennae  elements  N  is  obtained. 


Receiver  Performance 

Cost  function 


J(Y;  0)  =  log  det  F(Y\  0) 

where  F  is  the  Fisher  information  matrix 


fd2e(Y-,6)\ 

\  w  / 


which  is  strictly  positive  definite.  This  quantity  is  a  function  of 
the  quantities  rn.m,  and  thus  of  the  state  variables  r  and  p,  as  well 
as  the  signal  waveforms,  and  thus  of  the  code  Yn  m. 


-p.  l/i 


Performance  Optimisation 

The  cost  function  J  is  chosen  because  its  maximisation  implies 
minimisation  of  the  Cramer-Rao  lower  bound  on  the  parameter 
estimation  uncertainty.  There  are  four  constraints  to  be 
met.Firstly,  the  beamformer  weight  magnitudes  specify  the 
transmit  powers  for  each  antenna,  ie  \  Wn\2  =  pn  =  Qn,  with 
(li  =  1-  We  also  have  the  code  constraints 

■f'k  $n,k 

xj:  Tt  xn  —  0,  V  n  ^  k 
xj?  Ttt  xn  =  0,  V  n  7^  k  , 

where  T t  and  T tt  are  diagonal  matrices  specihed  by  the 
signalling  pulse,  and  the  xn  are  the  rows  of  X  (ie  the  sequence 
transmitted  on  antenna  n ).  p 


Conclusion 

space-time  coding  leads  to  a  new  design  concept  for  adaptive  radar 

estimates  are  unbiased  so  that  CRLB  can  be  utilised  for  performance 
measure 

processing  can  be  done  at  individual  antenna  elements  -  reduces  load 
on  central  processor  for  real-time  implmentation 

this  formulation  can  be  generalised  for  netcentric  radars 
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Open  Problems 

•  Are  the  constraints  compatible  ? 

•  How  do  we  find  the  best  code  ? 

•  Is  the  receiver  efficient?  -  Kershaw  and  Evans  proves  this 
for  SISO. 

•  How  does  one  maximize  //  in  the  presence  of  local  maxima? 


